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Abstract. A model of interstitial impurity migration is proposed which explains the redistribution 
of ion-implanted boron in low-temperature annealing of nonamorphized silicon layers. It is supposed 
that nonequilibrium boron interstitials are generated either in the course of ion implantation or at the 
initial stage of thermal treatment and that they migrate inward and to the surface of a semiconductor 
in the basic stage of annealing. It is shown that the form of the “tail” in the boron prohle with the 
logarithmic concentration axis changes from a straight line if the average lifetime of impurity interstitials 
is substantially shorter than the annealing duration to that bending upwards for increasing lifetime. 

The calculated impurity concentration prohles are in excellent agreement with the experimental data 
describing the redistribution of implanted boron for low-temperature annealing at 750 Celsius degrees for 
1 h and at 800 Celsius degrees for 35 min. Simultaneously, the experimental phenomenon of incomplete 
electrical activation of boron atoms in the “tail” region is naturally explained. 

1 Introduction 

The inhuence of the impurity diffusion mechanism on the form of dopant distribution after annealing of 5- 
doped layers has been investigated in m- The authors obtained a one-dimensional analytical solution of 
the system of equations that includes the conservation law for the immobile component (substitutionally 
dissolved impurity atoms) and the diffusion equation for the mobile impurity component [1]. It was 
supposed that exchange of impurity atoms occurs between mobile and immobile species. It was also 
supposed that the initial impurity distribution is described by the (5-function and that the generation 
rate of mobile impurity interstitials is proportional to the concentration of substitutionally dissolved 
impurity atoms. It follows from the solution obtained that the extended “tail”, formed in the region of 
low impurity concentration of dopant prohle after annealing, represents a straight line for the logarithmic 
concentration axis if the impurity atom migrates once. On the other hand, in the case of multiple events of 
migration of the interstitial impurity atom, the Gaussian distribution describes the shape of the impurity 
prohle in the “tail” region. 

As follows from the experimental data, after annealing the impurity concentration prohles can be sim¬ 
ilar to the Gaussian distribution. On the other hand, many of the measured impurity prohles, especially 
after low-temperature treatment, are characterized by an extended “tail” which looks like a straight 
line if the concentration axis is logarithmic (see, for example, the review of the experimental data in 


1 


0 )- Therefore, it was proposed in mm that the microscopic diffusion mechanism of impurity atoms 
be determined based on the form of the dopant distribution after annealing. It is worth noting that in 
mm continuous generation of mobile species in the region of high impurity concentration was supposed 
in the case where an impurity atom migrates once. In mm analytical solutions were obtained for the 
initial impurity distributions that, in contrast to mm, represent an uniform doped layer or the Gaussian 
distribution. It was shown that the “tails” also represent a straight line just as for the initial distribution 
described by (5-function mm- 

At present, the method of ion implantation with low budget annealing such as, for example, high- 
temperature spice annealing or low-temperature long-term thermal treatment, is widely used for fabri¬ 
cating modern semiconductor devices and integrated microcircuits. During ion implantation, a fraction 
of the implanted ions can occupy an interstitial position and participate in the interstitial migration at 
subsequent annealing. It is expected that this phenomenon will be most pronounced for implantation of 
light ions, for example, boron ions, and also for implantation of the medium doses of heavier ions 
cm“^) when the maximal concentration of impurity atoms is below the solubility limit but higher than 
the intrinsic carrier concentration at the annealing temperature |^. Indeed, in these two cases one does 
not observe the amorphization of an implanted layer. As a result, annealing of defects generated by ion 
implantation occurs during the total thermal treatment, instead of the solid phase recrystallization of 
the created amorphous layer at the initial stage of the thermal processing. Indeed, it is well known that 
during solid phase epitaxial regrowth (SPER) there occurs annealing of the main part of defects created 
by ion implantation. In addition, almost all implanted impurity atoms go over into the substitutional 
position, even for the concentration above the solubility limit. It will be expected that in the absence of 
SPER the impurity atoms that remained in the interstitial position will migrate to a significant distance 
before they become substitutionally dissolved. It is worth noting that in the case of low-temperature 
treatment with duration of ten minutes and longer, a significant number of the impurity interstitials can 
also be generated at the initial stage of thermal processing due to annealing or rearrangement of different 
kind of defects incorporating impurity atoms. 

Therefore, in contrast to the case of continuous generation of impurity interstitials during annealing 
described in mmmm, the purpose of this work is to investigate the characteristic features of ion- 
implanted impurity redistribution that can appear due to diffusion of the previously formed interstitial 
impurity atoms. 

2 The system of the equations describing interstitial diffusion 

Let us assume that the substitutionally dissolved impurity atoms, as well as the impurity atoms incorpo¬ 
rated into clusters, precipitates, and radiation defects are immobile and do not participate in diffusion due 
to the low annealing temperature. As the formation of extended “tails” due to the migration of impurity 
interstitials occurs in the region of low impurity concentration C < rij, we neglect the influence of the 
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built-in electric field on the dilfusion process. Here C is the concentration of substitutionally dissolved 
impurity atoms and rii is the intrinsic carrier concentration. Then, the system of equations that was 
proposed in mM and that describes migration of nonequilibrium impurity interstitial atoms and their 
subsequent transition to the substitutional position can be presented in the form given below. 

1. The conservation law for substitutionally dissolved impurity atoms: 

dC{x,t) _C^^{x,t) 

dt ^ ’ 

2. The equation of diffusion for nonequilibrium impurity interstitials: 
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where is the concentration of nonequilibrium interstitial impurity atoms (HA); and are 
the diffusivity and average lifetime of these impurity interstitials, respectively; is the generation 

rate of interstitial impurity atoms per unit volume due to the annealing of implantation defects and 
rearrangement of clusters (precipitates). It is worth noting that the total concentration of impurity atoms 
includes concentrations of the substitutionally dissolved impurity atoms, impurity atoms incorporated 
into clusters, precipitates, and radiation defects, as well as the concentration of impurity interstitial atoms 

G^A 


It is well known that interstitial impurity atoms can become substitutionally dissolved due to the 
recombination with vacancies (the so-called the Frank-Turnbull diffusion mechanism [^). In addition, 
a migrating interstitial impurity atom can kick-out the host atom and occupy a substitutional position 
(the reverse Watkins effect). Then, the average lifetime of IIA can be presented in the following form: 


^ (3) 

where k^^^ is the effective coefficient of recombination of interstitial impurity atoms with vacancies; 
is the vacancy concentration; k^^ is the coefficient describing transition of the interstitial impurity atom 
to the substitutional position due to the kick-out mechanism. 

It is worth noting that this paper investigates an extreme case of the redistribution of implanted atoms 
occupying interstitial position immediately after implantation. Therefore, the generation term G^^^ = 0. 
Then, Eqs. ([T]) and ([2]) are independent if the vacancy distribution is uniform {C^ {x,t) = const). If d^^ 
is also constant, one can obtain an analytical solution of Eq. ©• 


3 Analytical solution of nonstationary diffusion equation for in¬ 
terstitial impurity atoms 

An analytical solution of Eq. ([2]) can be obtained within the framework of the Green function approach 
[iniiniin]. For this purpose let us present Eq. ([^I in the following form: 
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where Iai = Vis the average migration length of impurity interstitials; 9 = tj is the relative 
time. 

Let us also assume that the distribution of nonequilibrium impurity interstitials after implantation 
and after rapid annealing of implantation defects incorporating impurity atoms is proportional to the 
distribution of the implanted atoms which can be approximately described by the Gaussian distribution. 
Then, the initial condition for impurity interstitials can be presented in the following form: 


C^\x) = C^\x, = 


-iAI 


AI 
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where is the maximal concentration of interstitial impurity atoms after implantation or at the initial 
stage of annealing; Rp and ARp are the average projective range of implanted ions and the straggling of 
the projective range, respectively. 

Let us obtain an analytical solution of the diffusion equation for nonequilibrium impurity intersti¬ 
tials (H]) in one-dimensional (ID) semiinfinite domain [0, -foo]. For the sake of simplicity, the Dirichlet 
boundary condition can be imposed on the surface and in the bulk of a semiconductor: 


^^^(0, 9) = ci 


AI 


C^^{+oo, 0 ) = 0 , 


( 6 ) 


where is the concentration of interstitial impurity atoms on the surface of a semiconductor. 

For this purpose, let us consider Eq. o in the infinite domain [—oo, -too] and use separation of 
variables [TOlIIl] assuming that the required solution is a product of two functions, one of which depends 
only on the relative time 9, whereas the other is the function of only the spatial coordinate x, i.e., 
C^^{x,0) = T{9)X{x). Then, = T'X, = TX" and Eq. (|4]) transforms into the equation 


T'X = l\jTX" - TX. 


( 7 ) 


Dividing the right- and left-hand sides of Eq. m by {l\jTX), we obtain the following equation: 

J2 rp Y P 

IaiI a. 

where A*^ is the parameter of division that depends neither on x, nor on 9, i. e., represents a certain 
constant. Using this property of A*^, one can obtain the following system of independent equations: 


T' 

^ _ ]2 \*2 


( 9 ) 


X" + A^A = 0, 


( 10 ) 


where A = A*^ — ■ 
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The solution of Eq. (jHl) has the following form; 


T{e)=Toexp[-l\jX*^0], ( 11 ) 

where Tq is the value of the function T{9) for the initial moment of annealing. 

The solution of homogeneous equation m can be obtained by employing the standard substitution 
X{x) = AQe°‘^, where Aq and a are some constants. Based on the obtained solutions of Eqs. (jH]) and 
and using the procedure described in m, it is possible to find the solution of the boundary value 
problem (jH), ([5]), and ([6]) for zero boundary condition on the surface of a semiconductor, = 0. This 
solution has the following form: 


POO 
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where G{x, 9) is the Green function 


G{x,^,9) = 
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To obtain the required distribution of the concentration of impurity interstitials for nonzero condition 
on the surface of a semiconductor, it is necessary to combine expression m with the well-known solution 
of the boundary value problem (|3]), ([S]) with zero initial condition ([S]): 


G^^{x, 9) = G^esi^, 9) + G^ 


■AI 
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A similar procedure can be used to obtain a solution for the reflecting boundary condition on the 
surface and zero Dirichlet boundary condition in the bulk of a semiconductor: 




dx 


= 0 , 


G^\+oo, 9) = 0. 


x—0 


This solution has the following form: 


pOO 

G^\x,9)= G{x,^,9)G^\Odi, 

JO 


(15) 
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where the Green function is 


G{x,^,9) = 


exp 


Al\,9 


exp 
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20^ 

The explicit expressions for the integrals in m and m can be obtained using the modern software 
allowing symbolic computation, such as “Mathematica” El, “Mathcad” m, and others. Unfortunately, 
the expressions obtained are too lengthy to be presented in this paper. 

The total concentration of impurity atoms after annealing can be calculated by Eq. ©■ Integrating 
this equation, we obtain 
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C{x, t) = f C^\x, t) dt + Co{x), 
Jo 


( 18 ) 


where Cq (x) is the initial distribution of substitutionally dissolved impurity atoms at the initial stage of 
annealing that can be approximated by the Gaussian distribution: 


Co{x) = Cm exp 


{x - Rp)' 

2ARI 


(19) 


Here Cm is the maximal concentration of substitutionally dissolved impurity. 

Unfortunately, within the framework of the software for the above-mentioned symbolic computation 
it is not possible to obtain an explicit expression for integral (fTSl) with C^^{x,t) described by expression 
da or (HU). Therefore, the calculation of integral (jl8l) was carried out by an approximate numerical 
method using the quadrature formula (Gauss 16-point formula) [161117) . 


4 Simulation of the interstitial redistribution of boron implanted 
into silicon during low-temperature thermal treatments 

In Fig. [T] the calculated profiles of ion-implanted boron after low-temperature annealing are presented 
which show the change in the shape of the low concentration “tail” formed due to the interstitial diffusion 
for different ratios between the average lifetime and the duration of thermal treatment tann- It is 
supposed that silicon is implanted with 1 keV boron ions to a dose of 1.8x10^^ ion/cm^ {Rp = 0.0056 
/im; Ai?p = 0.0038 /rm eg) and then subjected to annealing with duration of 10 min at a temperature 
of 750 °C. The average migration length of nonequilibrium boron interstitials has been chosen to be 
equal to 0.08 ^m. It is clearly seen from Fig. [T] that if the average lifetime of impurity interstitials 
is significantly shorter than the annealing duration tann, the “tail” in the boron concentration profile 
represents a straight line in the case of the logarithmic axis of concentration. With increase in the average 
lifetime the shape of the “tail” becomes convex upwards, being characterized by an increasing slope 
angle in the bulk of a semiconductor similar to the Gaussian distribution. Thus, the calculations presented 
in Fig. [T] show clearly that for nonzero initial distribution of interstitial impurity atoms the “tail” in the 
low concentration region of impurity distribution has the form of the Gaussian distribution even for one 
event of impurity interstitial migration if > tann, whereas in the case of continuous generation of 
impurity interstitials it occurs only in the case of numerous events of interstitial migration [Dig¬ 
it follows from the experimental data that the “tail” in the boron concentration profile after low- 
temperature treatment of ion-implanted layer can have a shape of a straight line muni HD as well 
as a convex form characteristic for the Gaussian distribution [i m m [25] if the concentration axis is 
logarithmic. For example, in Fig. [2] the boron concentration profile is presented which was measured in 
[6] by a method of the secondary ion mass spectroscopy (SIMS). Both the p-type and n-type Czochralski 
grown (CZ) commercial silicon wafers with conductivities of 4-6 Hem were used. Boron was implanted at 
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Figure 1: Calculated boron concentration profiles after low-temperature thermal treatment of implanted 
silicon substrate for different ratios between the average lifetime and the duration of annealing tann 

room temperature with an energy of 40 keV to a dose Q of 3x 10^'*’ cm“^. SIMS profile of the implanted 
boron after annealing at 750 °C with duration tann = 1 h is presented in Fig. [5J by filled circles. It 
can be seen from Fig. [2]that the maximal boron concentration after implantation Cm = 3.04x10^ ^m“^ 
is near the limit of boron solubility in silicon Csoi = 2.33x10^ /im“^ for a temperature of 750 °C [26] . 
According to m. the boron diffusivity of substitutionally dissolved atoms Di for this temperature is 
equal to 2.658x10“^° /rm^/s. Then, for the thermal treatment for 1 h the characteristic diffusion length 
of impurity redistribution is Li = \/Ditann = 9.78x10“^ /rm or it is approximately equal to 1 nm. On 
the other hand, it follows from Fig. |2|that the actual Li is approximately equal to 0.1 /rm, i.e., 100 times 
longer. If we suppose that boron diffusion occurs due to the formation, migration, and dissociation of the 
“impurity atom - intrinsic point defect” pairs, it is necessary to accept that the average concentration of 
nonequilibrium point defects responsible for impurity diffusion is approximately 10'^ times higher than the 
thermally equilibrium value. Not rejecting the possibility of such a strong radiation-enhanced diffusion, we 
shall carry out modeling of the redistribution of ion-implanted boron based on a more natural mechanism 
of migration of nonequilibrium impurity interstitials. 

The following values of the model parameters were used to provide the best fit of the calculated boron 
concentration profile to the experimental one: The parameters prescribing the initial distribution 
of implanted boron: Q = 3.2x10^^ cm“^; Rp = 0.145 /im; ARp = 0.042 /rm; the fraction of the 
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Figure 2: Calculated boron concentration profile (solid line) after thermal treatment of implanted silicon 
substrate at 750 °C for 60 min. Gaussian distribution (dotted line - total boron concentration and 
dashed line - concentration of interstitial boron atoms) was used to approximate the initial impurity 
profile. The dash-dotted curve represents interstitial boron after annealing. The experimental data 
(impurity distribution after annealing represented by filled circles) are taken from Ref. [B] 

boron atoms in the interstitial position = 34.2 %. The parameters specifying the process of 
interstitial diffusion: the average migration length of boron interstitials Iai = 0.11 pm. The reflecting 
boundary condition is imposed on the surface of a semiconductor when describing diffusion of boron 
interstitials. It is also supposed that the average lifetime of nonequilibrium impurity interstitials is 
equal to the duration of annealing tann = 1 h. It is worth noting that the parameters prescribing the 
initial distribution of implanted boron is approximately equal to the values tabulated in [18j : Rp = 0.131 
pm; ARp = 0.0451 ^m; Sk = - 0.8 pm; Rm = 0.143 pm, where Sk and Rm are respectively the skewness 
and the position of a maximum of impurity distribution as implanted which is described by the Pearson 
type IV distribution [18] . 

It can be seen from Fig. [2| that there is an excellent agreement between the calculated total boron 
concentration profile and boron profile measured by SIMS. A little difference in the near surface region 
follows from the assumption about the reflecting boundary for boron interstitials. In point of fact part 
of impurity interstitials cross over the surface of a semiconductor. 

The modeling results of interstitial diffusion of ion-implanted boron during a low-temperature anneal¬ 
ing at a temperature of 800 °C for 35 min are presented in Fig. [3| The experimental data of [23] were 







used for comparison. In the Ref. [53] boron was implanted with an energy of 70 keV to a dose Q of 10^® 
cm“^. For the energy used in ion implantation, the parameters of the Pearson type IV distribution are: 
Rp = 0.219 /xm; ARp = 0.0606 ^m; Sk = -0.9 [18]. It is worth noting that for asymmetry Sk = -0.9 
the position of a maximum of boron concentration R^ = 0.236 fim is close to the value Rp = 0.219 fim. 
It means that, just as in the previous case, one can neglect the asymmetry of impurity distribution and 
describe the boron profile after implantation, including the distribution of boron atoms in the interstitial 
position, by the symmetric Gaussian distribution ([5]) |28j . 

The diffusivity of the substitutionally dissolved boron atoms Di is equal to 1.928x10“® /xm^/s for 
a temperature of 800 °C m- Then, the characteristic diffusion length of impurity redistribution Li = 
V Ditann = 2.01x10“® /im Or approximately 2 nm for thermal treatment during 35 min. On the other 
hand, it can be seen from Fig. [3] that the actual value of Li is approximately 0.1 /xm, i.e., 50 times 
longer. To provide this value of L^, the time-average concentration of the nonequilibrium point defects 
responsible for the impurity diffusion should be 2.5x 10® times higher than a thermally equilibrium value. 
It is necessary to note that in the case of radiation-enhanced diffusion of impurity atoms, due to the 
formation, migration, and dissociation of “impurity atom - intrinsic point defect” pairs, all impurity 
atoms in the “tail” region should be substitutionally dissolved, i.e., to be electrically active. However, 
the experimental data of [23] show that the significant fraction of boron atoms is electrically inactive in the 
“tail”. This phenomenon is easily explained within the framework of the interstitial diffusion mechanism 
described above if the lifetime of the interstitial boron atoms is near or longer the annealing duration 
tann- Indeed, it follows from the later condition that a significant part of boron atoms are staying in the 
interstitial position, i.e. they are electrically neutral. 

The following values of the model parameters were used to provide the best fit of the calculated 
boron concentration profile to the experimental one (see Fig. |3]): The parameters prescribing the 
initial distribution of implanted boron: Q = 1.12x10^® cm“^; Rp = 0.24 /xm; ARp = 0.054 /xm; 
the fraction of the boron atoms in the interstitial position = 26.5 %. It can be seen in Fig. [3|that 
the maximal boron concentration after implantation Cm = 8.3x10^ /xm“® is more than 2 times higher 
than the limit of boron solubility in silicon Csoi = 3.431x10^ /xm ® for a temperature of 800 °C [26] . 
The parameters specifying the process of interstitial diffusion: the average migration length of 
boron interstitials Iai = 0.14 /xm. As in the previous simulation, the reflecting boundary condition is 
imposed on the surface of a semiconductor when describing diffusion of boron interstitials. The average 
lifetime of nonequilibrium impurity interstitials = 48,3 min is chosen to be longer than the duration 
of annealing tann = 35 min. It is worth noting that the parameters prescribing the initial distribution of 
implanted boron is approximately equal to the values listed in the above-mentioned tables [18] . 

It can be seen from Fig. [3] that there is excellent agreement between the calculated total boron 
concentration profile and boron profile measured by SIMS. It is interesting to note that the concentration 
of interstitial boron atoms is approximately equal to the total boron concentration at the end of the “tail”. 
It means that this part of the “tail” is electrically inactive that is also consistent with the experimental 
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Figure 3: Calculated boron concentration profile (solid line) after thermal treatment of implanted silicon 
substrate at 800 °C for 35 min. Gaussian distribution (dotted line - total boron concentration and 
dashed line - concentration of interstitial boron atoms) was used to approximate the initial impurity 
proHle. The dash-dotted curve represents interstitial boron after annealing. The experimental data 
(impurity distribution after annealing represented by filled circles) are taken from Ref. [23] 

data [53]. 

This paper considers the diffusion of boron interstitial atoms generated either during ion implantation 
or at the initial stage of annealing. In this meaning the case of interstitial diffusion under consideration 
is opposite to the investigations of dJIlllIS] where continuous generation of impurity interstitials is 
supposed. It follows from the simulation results obtained (see Figs. 2 and 3) that the model of diffusion 
of previously generated boron interstitials allows one to explain the formation of “tails” in nonamorphized 
silicon layers during low-temperature annealing with duration of 1 h or shorter. On the other hand, some 
increase in the boron concentration in the “tail” region occurs for annealing duration of 21 h at 800 °C 
[23] in comparison with the simulated boron profile for 1 h annealing presented in Fig.[3j We suppose that 
the additional generation of boron interstitials for such long-time treatments occurs due to the annealing 
or rearrangement of clusters or radiation defects. 
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5 Conclusions 


To explain the redistribution of ion implanted boron during low-temperature annealing of nonamorphized 
silicon layers, a model of interstitial impurity migration is proposed. In contrast to other models of in¬ 
terstitial diffusion, it is supposed that nonequilibrium boron interstitials are generated either during ion 
implantation or at the initial stage of annealing and migration inward and to the surface of a semicon¬ 
ductor during the basic stage of annealing. It was shown that such kind of interstitial diffusion results 
in a “tail” in the boron concentration profile which has a form of a straight line for the logarithmic con¬ 
centration axis if the average lifetime of impurity interstitials is significantly shorter than annealing 
duration tann and becomes convex upwards similar to the Gaussian distribution with increasing . 

The calculated impurity concentration profiles are in excellent agreement with the experimental data 
that describe redistribution of implanted boron for low-temperature annealing at 750 °C for 1 h and at 
800 °C for 35 min if is equal to 1 h and to 48.3 min, respectively, i.e., the average lifetime of boron 
interstitials is equal or rather longer than the duration of annealing. The average migration length of 
impurity interstitial atoms is equal to 0.11 /rm and to 0.14 /rm for the annealing at temperatures 750 and 
800 °C, respectively. The simulation results also show that 34.2 % and 26.5 % of the implanted boron 
atoms occupy interstitial position at the initial stage of annealing. 

It is worth noting that the experimental phenomenon of incomplete electrical activation of boron 
atoms in the “tail” region is naturally explained within the framework of the model proposed. This 
fact and the excellent agreement with the measured boron proHles give evidence in favor of the model 
considering the formation of the main part of boron interstitials generated either during ion implantation 
or at the initial stage of annealing if we deal with low-temperature treatment of nonamorphized silicon 
layers. 
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